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ABSTRACT: Partially quenched theories are theories in which the valence- and sea- 

quark masses are different. In this paper we calculate the nonanalytic one- loop corrections 
!>■ ■ 

Q\ • of some physical quantities: the chiral condensate, weak decay constants, Goldstone boson 



masses, Bj{ and the K + — ► 7r + 7r^ decay amplitude, using partially quenched chiral pertur- 
bation theory. Our results for weak decay constants and masses agree with, and generalize, 
Ph, 

results of previous work by Sharpe. We compare Br and the K + decay amplitude with 

i-G ■ 

their real-world values in some examples. For the latter quantity, two other systematic ef- 
fects that plague lattice computations, namely, finite-volume effects and unphysical values 
of the quark masses and pion external momenta are also considered. We find that typical 
one-loop corrections can be substantial. 
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1. Introduction 



Recently, large scale numerical lattice QCD computations have started to move away 
from the quenched approximation by taking virtual-quark loop effects into account. However, 
since the computation of quark determinants is much more expensive than that of quark 
propagators, one is often restricted to only very few values of the sea-quark mass, while, 
given an ensemble of "dynamical fermion" gauge configurations, one can explore a much 
larger range of values for the valence-quark masses. This naturally leads us to consider 
partially quenched theories, which are theories in which the valence- and sea-quark masses 
are not equal. 

Chiral Perturbation Theory (ChPT) plays an important role in the analysis of results 
from lattice computations (see for instance ref. [1]). Since the quenched or partially quenched 
approximations of lattice QCD are theories different from the full unquenched theory, ChPT 
needs to be adapted to these situations. For the quenched case, this has been extensively 
investigated (for a rather complete list of references, see ref. [2]), but much less work has 
been done for the partially quenched case. 

Partially Quenched Chiral Perturbation Theory (PQChPT) was developed in ref. [3], 
and used to calculate one-loop expressions for Goldstone boson masses and decay constants 
in ref. [4]. In this paper, we revisit Goldstone boson masses and decay constants, and we 
also give the one-loop expressions for the chiral condensate, and the K + — > 7r + 7r° decay 
amplitude. (For an application to heavy-light decay constants and B-parameters, see ref. 
[5]-) 

In the partially quenched case (unlike in the quenched case), the rj' meson is heavy in the 
sense that its mass does not vanish in the chiral limit. One can therefore approach PQChPT 
in the same way as one does in the unquenched case, where the rj is integrated out, and 
only the Goldstone bosons are kept in the effective lagrangian. This is the approach taken 
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in ref. [4]. In the completely quenched case, however, the rj' is light, and needs to be kept 
in the effective lagrangian, along with the other Goldstone bosons. The theory depends on 
a mass scale mo, independent of the quark masses, which in the unquenched theory is the 
singlet part of the rf mass, but in the quenched theory no longer appears in any pole mass 
[6,7,8,9]. A systematic expansion in the quenched case only exists if we treat the ratio of 
mo to the chiral-symmetry breaking scale as an independent small parameter, in addition to 
the light quark masses, and if we stay away from the chiral limit, where infrared divergences 
occur [8,9,10,11]. 

In partially quenched lattice computations, it may happen that the scale set by the 
sea-quark mass is not small compared to the scale mg. In this paper we therefore keep the 
rf in the effective theory, and show that one can still develop PQChPT systematically if 
one assumes, as in QChPT, that the ratio of mg to the chiral-symmetry breaking scale is 
reasonably small. This generalizes the results of ref. [4]. 

The outline of this paper is as follows. In section 2, we give a quick review of PQChPT, 
and discuss the role of the rf in more detail. We explain more precisely how our calculation 
generalizes that of ref. [4]. We then calculate in section 3 the chiral condensate and Gold- 
stone boson masses and decay constants to one loop as a function of the quark masses, and 
investigate their dependence on ttiq. We restrict ourselves to the theory with N degenerate 
sea quarks. In section 4, we calculate Br for non-degenerate valence-quark masses and the 
K+ — > 7r + 7r^ decay amplitude for degenerate valence-quark masses. For K + — > n^ii® we also 
discuss other systematic errors which have affected lattice computations of this quantity to 
date. In section 5, we give numerical examples of the role of one-loop corrections for and 
K+ — > 7T + 7r° in a comparison between the real world and the partially quenched theory with 
parameters typical of a lattice computation. More details on the assumptions underlying 
these numerical examples will be given in due course. We end with our conclusions. 
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2. Essentials of Partially Quenched Chiral Perturbation Theory 



Consider a QCD-like theory with n flavors of quarks q^, each with arbitrary mass m^. 
We then partially quench the theory by adding n — N flavors of (unphysical) quarks 
obeying bosonic statistics [12], which we will call ghost-quarks. The ghost-quark masses are 
chosen to equal those of the first n — N physical quarks. These n — N quarks are then 
quenched since their loop contributions are exactly canceled by their bosonic counterparts. 
The remaining iV quarks contribute through (sea-) quark loops to physical quantities. With 
the n — N quarks identified as valence quarks, the theory corresponds to partially quenched 
QCD with N flavors of sea quarks. Recent partially quenched lattice computations use 
degenerate sea-quark masses, hence, in the following, we will set all sea-quark masses raj, 
i — n — N + 1, • • • , n, equal to mg. 

The full chiral symmetry of the theory is the semi-direct product of graded groups 
G = [SU(n\n — N) <g> SU(n\n — N)\ (s)U(l) after the anomaly has been taken in account 
[8]. We will briefly discuss the construction of the lagrangian for PQChPT for our choice of 
quark masses. 

The unitary field E is defined through 

E = exp(2^//) (1) 
from the (2n — N) x (2n — N) hermitian matrix field 

where is the n x n matrix of ordinary mesons made from the n ordinary quarks and 
antiquarks, <f> is the corresponding (n — N) x (n — N) matrix for ghost-quark mesons, and 
X is an (n — N) x n matrix of fermionic mesons made from a ghost-quark and an ordinary 
antiquark. / is the tree-level pion decay constant. The (2n—N) x (2n—N) quark-mass matrix 

4 



M is defined by = m^Sij, where, as already discussed, the masses for i — 1, • • • , n — N 
are arbitrary, and equal to the ghost-quark masses, i — n + 1, • • • , 2n — N . The remaining 
N masses i — n — N + 1, • • • ,n, are degenerate and equal to mg. 

As discussed in refs. [8,3], the super- rf field $o = str($) is invariant under the full chiral 
group G and introduces new parameters 7Bq and a into the O(jo 2 ) euclidean lagrangian for 

PQChPT, which reads 

£ 2 f 2 2 

C = L str (^S^st) - LE str (mS + M S t) + ^ $g + | (^$ ) (d M <t> ) . (3) 



We define 



M?j = /iK + mj) (4) 



(for % — j — n — N + 1, • • • , n, this simplifies to = 2fimg = Mg^). It is instructive to 

display the two-point functions for the neutral mesons (i.e. the diagonal fields explicitly. 

In the diagonal basis with states « = 1, • • • , 2n — N, corresponding to q\q\, <72<?2> ' ' ' > 

and their ghost-quark counterparts, these two-point functions are, in momentum space [3], 

r M V» 1 (ml + ap 2 )(p 2 + Ml s ) 

" ? + M l ( 3 + iV «) (P 2 + M2 )(P 2 + M2 ,)(P 2 + ™ 2 ') ' 



where 



and 



+ 1, for 1 < i < n 

-1, for n + 1 < i < 2n -N ' 



(6) 



2 _ Mfg + iVmjj/3 
"V " 1 + iV«/3 ' (7) 

which is the square of the pole mass for the $o two-point function, as can easily be veri- 
fied from Eq. (5). It is equal to the square of the "77'" meson mass in the SU(N) theory 
constructed only from the unquenched quarks [3]. It does not vanish in the chiral limit 
mg — > 0. 

A simplification occurs in partially quenched QCD. For a correlation function in which 
only k out of the n — N valence quarks are on the external lines, the rest of the n — N — k 
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valence quarks do not contribute at all, and therefore this correlation function does not 
depend on their masses. In most of this paper, we will only consider quantities involving at 
most two valence quarks, i.e. we will only consider valence quarks qi with i — 1,2. 

Compared to quenched ChPT (QChPT), PQChPT introduces an extra mass scale Mgg. 
We can therefore consider various possible "expansion schemes." Let A p be the cutoff of the 
partially quenched theory. The first scheme, adopted in ref. [4], takes M\\, M22, Mgg < 



rrirf ~ 



Ap. This corresponds to the usual set-up for power counting in ChPT, in which the 77 



is not a Goldstone boson but instead a "heavy" particle, due to the presence of the scale mo, 
which does not vanish in the chiral limit. In this scheme, the field $0 can be integrated out, 
absorbing dependence on the parameters a and itlq into 0(p 4 ) terms in the chiral lagrangian. 
This can be illustrated by considering for instance the two-point function G\^\ (cf. Eq. (5)). 
We first write the second term in G\ \ as 



i 2 + ap 2 ( ml,-M 2 St 



3 + Na V( m J/- M ll) 2 



1 1 



p 2 + M\ x p 2 + m 2 , 



Mg S — M 2 X 1 
+ m 2 , - M\ x (p 2 + M 2 n ) 2 1 i;<) 



For large m^, the expansion of the term containing the pole l/(p 2 + m 2 ,) in powers of 
M 2 ^/m 2 , and Mgg/m 2 ; generates contact terms which can be absorbed into 0(p 4 ) coeffi- 
cients of the chiral lagrangian. Similarly expanding the remaining terms in Eq. (8) gives the 
leading terms 

+ ■ (9) 



The poles in this expression, which no longer depends on a and to 2 ,, lead to nonanalytic 
terms at one loop. For details and further discussion, see ref. [4]. 

In QChPT, there is no heavy particle with a mass determined by the parameter 7Bq, 
as can been seen by setting iV = in Eq. (5). (The $q two-point function vanishes for 
N = [8].) The second term in Eq. (5) is proportional to + ap 2 , and has a double 
pole at the Goldstone meson masses M 2 ^. (Of course, no Mgg dependence remains.) Hence, 
this term needs to be kept in the quenched approximation, and the decoupling as described 
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above does not occur. It turns out that, in order to make sense of the chiral expansion, 
we need to assume that the parameter (mQ/3)/(47r/ 7r ) 2 is small, in addition to the usual 
requirement that M^/(47r/ 7r ) 2 be small [8,9]. QChPT therefore corresponds to a second 
"scheme" in which M\\, M22, mo < h-p and Mgg > A p which corresponds to freezing out 
sea-quark loop effects, effectively setting N = 0. We note that in the unquenched theory 
(m§/3)/(47r/ 7r ) 2 w 0.09 (where f n = 132 MeV is the physical pion decay constant). It is 
believed that the value of mg is not very different in the quenched theory [1]. 

Since the partially quenched theory "interpolates" between the quenched and unquenched 
theories, it is natural to consider a third scheme, which we will now explain. First, as in ref. 
[4], we take M\\, M22, Mgg < A p , and M\\, M22 < m^. However, we would like to leave 
the ratio of Mgg and mg arbitrary. Let us consider the possibilities. If Mgg < m^i we can 
systematically study the effective theory for mesons with masses M\\, M22 and Mgg by in- 
tegrating out the rj', because now M\\, M22, Mgg < m^, as in ref. [4]. If, on the contrary, 
Mgg is of order uiq or larger, this would be incorrect. Now, we can only systematically 
investigate the effective theory for mesons with masses Mjj, M22 and Mgg, if we keep the 
7/ in the effective theory, and, as in the completely quenched case, assume that mp < A p . 

In partially quenched lattice computations, often Mgg is fairly large, and we may well 
be in the situation that indeed Mgg is of order thq. Therefore, in this paper, we will be 
interested in calculating the nonanalytic dependence on the quark masses M\\, M22, Mgg 
for various quantities at one loop, taking Mjj, M22, Mgg, uiq < A p and M\\, M22 < m^i. 
Note that we can still consider the case that Mgg < mg if M\\, M22 < moi which leads 
us back to the assumptions made in ref. [4]. Indeed, when we expand our results for meson 
masses and decay constants in Mgg/m^, we obtain those of ref. [4].* On the other hand, 
if we take Mgg large and expand in Mn/Mgg, M22/Mgg and mo/Mgg, we obtain the 

* Due to a difference in normalization, our tree-level weak decay constant / is by a factor \/2 
larger than that of refs. [4,9]. This should be taken into account before comparing results. 
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quenched (i.e. N = 0) results of refs. [8,9,13]. 

Before we present our results, we should address one more issue. In our scheme, there 
will be one-loop contributions proportional to logm^, coming from $q tadpoles. For values 
of the meson masses which are small relative to m^/, such terms can be absorbed into the 
0(p 4 ) coefficients (after expanding in M^/m^i). However, this is not true for Mgg when 
Mgg is of order mg or larger. In this case, the dependence of $o tadpoles on Mgg will be 
complicated. Moreover, the contributions from $q tadpoles will depend on other r\ coupling 
constants (analogous to TOq and a) which are basically unknown. (We did not include these 
couplings in Eq. (3); for a more complete expression, see refs. [8,3].) 

Therefore, in this paper, we will take the point of view that we keep Mgg fixed, and 
calculate the nonanalytic dependence on the valence-quark masses. This, then, allows us to 
ignore contributions from <3>o tadpoles, which we will do in the rest of this paper. It does 
not affect the coefficients of the other chiral logarithms. Note that we do not assume the 
ratios of the valence-quark masses and the sea-quark mass to be small; our scheme allows 
for arbitrary values of these ratios. 



3. Chiral condensate, masses and weak decay constants 



We list the one-loop expressions for the condensate, meson masses and decay constants 
for a theory with nondegenerate valence-quark masses mi and 1122 and sea-quark mass mg 
in terms of the bare parameters /, M 2 -, itlq and a, including only nonanalytic terms (with 



M? s = (M? i + M 2 ss )/2,cf.E q . (4)): 



M 2 -f 2 

K<9*9*)]l-loop = l j— 



2N 



M 



iS 



(47T/)' 



log 



A2 



3(4tt/) ; 



M 2 - AMI + (m 2 - 2AMt) log 



a [-^mIi— loopt/wll— loop 5 



(10) 



M 2 2 



= M 2 9 I 1 + - 

1-loop 12 V 3(4tt/) 2 



2 



M ll O^ 2 - AM ll) M ll 
A*f 2 " M l 1 l0g ^T 

M| 2 (.M 2 - AMj 2 ) M| 2 
M| 2 - M 2 X 



;n) 



[/l2ll-Iopp // = 1 " ^2 ( M 1S lo S + M 2S j (12) 



p "-p 



where 



£s »g_^l = ^22 -11 (13) 



m 2 -mi _ M| 2 - 
m2 + mi M| 2 + M 2 : 

All dependence on itlq and a is embodied in the quantities M 2 and A, with 



u2 _ M SS m l _ 3y 2 

^ = M| 5 + Nm 2 /3 -TTi^ M ^' 



aM|^ + a + 3Ny< 

(M 2 ss + Nm 2 /3) 2 ~ (1 + Ny)'< 



A = — ggT— 0/ " = ^L^y , (is) 



where we introduced the ratio y = (m^/ByM 2 ^. The results of ref. [4] (with 0:4 = 05 = 
«6 = «8 — in Eqs. (13-20) of ref. [4]) can easily be obtained by expanding M 2 and 
A in Mgg/m 2 , (or mJ^/ttiq), and keeping the leading order terms, SMgg/N and 3/iV, 
respectively. (Subleading terms are of higher order in the chiral expansion, and can be 
dropped.) This corresponds mathematically to taking the limit y — > 00. Let us comment 
briefly on this comparison between our results and those of ref. [4]. In our case, we keep the 
7/, whereas in ref. [4] the 7/ is integrated out. One would therefore in general expect that 
in order to "match" the two theories, we would need to adjust the bare parameters. For the 
quantities considered here, it turns out that, at one loop, all nontrivial adjustments come 
from <3>o-tadpole contributions. However, we did not have to, and hence did not, include such 
contributions, as was explained in section 2. We conclude that, for the Goldstone meson 
masses and decay constants, no adjustment is needed. 
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The results of QChPT [8,9] can be obtained by taking N = 0, for which M = m 2 , and 
A — a. 



With degenerate quark masses my = mj = 777.2 m Eq. (11) and Eq. (12), we obtain 




As in QChPT, the ratio of Myy to its tree-level value, Myy j is singular, while fyy is 
regular in the chiral limit my — > 0. Moreover, fyy does not depend on 777q and a. 



We can define Mya and fyg to one loop from Eqs. (11,12) by replacing m\ by my 
and r?72 by mg. In the chiral limit my — > 0, keeping 7775- fixed, we have e — > 1 and 
the nonanalytic functions multiplying .M 2 and A in the expression for /yg, cf. Eq. (12), 
diverge. It is obvious from Eq. (11), that the ratio of Mya to its tree-level value is regular 
in the chiral limit, my — > 0. 

Of course, all quantities considered here also receive analytic contributions from 0(p 4 ) 
terms in the chiral lagrangian. Since they are obtained from tree-level diagrams with 0(p 4 ) 
vertices, they do not explicitly depend on the parameters 777q and a, and therefore are 
identical to those reported in ref. [4], to which we refer for their explicit form. In ref. [4] it 
was pointed out that the quantity 

Myy — Myig with my/ = 2my — mg , (18) 

is independent of 0(p 4 ) coefficients. 

Another quantity which is independent of 0(p 4 ) coefficients, first introduced in ref. [14] 
for the quenched case, and also considered in ref. [4] for the partially quenched case at 
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y — oo, is 



(19) 

For any rug, this quantity diverges when mi — » with m2 fixed (ie. for e — > 1). Note that 
[/l2]l-loop i s an even function of e (because of symmetry under the interchange mi <-> 7712), 
so that [ /12/V/11/22 ]l-loop = 1 + ^(e 2 ) for e -> 0. 

From Eqs. (11,12,16), we see that the coefficients of the chiral logarithms of the valence- 
quark mass depend on the ratio y of the parameter and the sea-quark mass, through the 
quantities M 2 and A (cf. Eqs. (14,15)). From (partially) quenched lattice data, it is esti- 
mated that presumably has a value m 2 K /2 ^ mg/3 ^ m 2 K (m^ = 496 MeV is the phys- 
ical kaon mass) [1,15]. Typical lattice computations have N = 2 and m 2 K ^ ^ 2m 2 K . 
These values of m^ and Mgg correspond to y ranging from y«l/4toy~l. This leads 
to M 2 /M^ s = 1/2 for y — 1/4 and to M 2 /M^ s — 1 for y — 1. For y — > 00 one obtains 
_A/f 2 /Mjg = 3/2. This shows that for relatively heavy sea quarks, there is a clear depen- 
dence of the coefficient of the chiral logarithms on mg. (Experience with quenched lattice 
data [1] indicates that it is hard to fit the chiral logarithms reliably, partially because of the 
"competition" of 0(p 4 ) coefficients. This may make it difficult to see the y dependence of 
the chiral logarithms in practice.) 

The quantity A also has an effect on the coefficients of the chiral logarithms, in particular 
for values of the valence-quark mass of order of the sea-quark mass. Taking again N = 2, we 
find A = (8a + 3)/18 for y = 1/4 and A = [a + 6)/9 for y = 1, while A = 3/2 for y -> 00. 
It is clear that A is more sensitive to the value of a for smaller values of y. We note that for 
my/mg = 1, M 2 ^/m 2 , = 1/(1 + Ny) (for a = 0), so that our results may not be reliable for 
smaller values of y. 
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4. Bx and K + — > 7r+7r° decay 



In this section, we will generalize earlier quenched one-loop calculations for [9,16,13] 
and K + — > 7r+7r° [13] to the partially quenched case. In the following, u, d and s denote 
valence quarks with masses mj, mi and rri2 respectively. The kaon B-parameter B% is 
defined as (with Mr the mass of the ds meson, i.e. at tree level Mr = M12) 

in which the four-quark operator is defined by 

(QiQjQkQ^LL = (QiLl^QjL) (QkLl^QlL) , ( 21 ) 

where = ^(1—75)5 is a left-handed quark field. The denominator in Eq. (20) is the matrix 
element {K \(sdsd)^^\K°) evaluated by vacuum saturation. The 0(p 2 ) weak-interaction 
operator in ChPT corresponding to (sdsd)n, which is the AS = 2 component of a 27-plet 
under SU (3)^, is 

O' = a 27 fg(E0"Et) /(E9 M Et)/ , (22) 
where the tensor t 1 ^ is defined by setting 

41 = 1, (23) 

while all other components are equal to zero. The parameter a 2 7 is the only new 0(p 2 )- 
operator coefficient; its value is determined by QCD dynamics. At tree level, one obtains 
[17] 

(K°\O'\K ) = , (24) 

and therefore, 

B K = ^ = B (25) 
12 



Wk = f a ^ tree level). Since the partially quenched theory is different from unquenched 
QCD, the partially quenched value of the coefficient a 2 7 is i n principle different from the 
QCD value. We will make this explicit by using a subscript or superscript p to denote bare 
parameters of the partially quenched theory, specifically, of 7 , B p and f p . 

The partially quenched one-loop result for B P K with non-degenerate quark masses, keep- 
ing only nonanalytic terms, is 

B p K = BPx (26) 



1 + ( - 2 ( 3 + e ) lQ g "aF - ( 2 + e ) l0 § C 1 - e ) - 3e 



2 

+ 



3(4tt/ p )< 



^(_log— - 

+ AMf 2 ( 2 + e 2 - 1 - 2e2 - £3 log 1±£ + 2 e 2 log ^ M - 



A 2 

where e is defined in Eq. (13). It can easily be seen that Eq. (26) does not depend on wlq and 
a in the case of degenerate quark masses (e = 0), just as in QChPT. Actually, apart from 
changing the values of B and / to their (partially) quenched values, (partial) quenching does 
not introduce any change in the nonanalytic one-loop corrections of Bj^ in the degenerate 
case. For the quenched case this was already discussed in ref. [9], and the fact that this is 
also true in PQChPT does not come as a surprise, since PQChPT is "in between" ChPT 
and QChPT (see also the discussion of K + — > 7r + 7r° below). In the nondegenerate case, the 
partially quenched result follows from the quenched result by replacing Wq — > M. , a — > A. 
For a discussion of contributions from 0(p 4 ) coefficients, see ref. [13]. 

The AS = 1, AI = 3/2 K + — > 7r + 7r° decay amplitude is proportional to the weak 
matrix element 

(tt+tt ! (sduu + suud - sddd) LL \K + ) . (27) 

The four-fermion operator is the AI = 3/2 component of the same 27-plet that also contains 
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the operator (sdsd) LL (Eq. (20)). To 0(p 2 ) in ChPT, the operator is represented by 

4 = a 27 rg(E^Et) i *(E^Et) i « , ( 28 ) 



where the tensor r. . has nonzero components 



21 12 12 21 1 



r 31 = r 13 = r 31 = r 13 
,22 22 ' 



2 ' 



'32 = rM = — (29) 
(all other components vanish). The parameter a 27 is the same as in Eq. (22). The aim is 
then to calculate the matrix element (7r + 7r |O4|K + ) to one loop. 

The lattice determination of (7r + 7r°|04|K + ) was compared with its real-world value in 
great detail at one loop in ChPT and QChPT in ref. [13], and here we will only discuss what 
is new in the partially quenched case. 

All attempts to compute this matrix element on the lattice have been restricted to 
the mass-degenerate, quenched theory, and moreover, all mesons are taken to be at rest 
[18,19,20]. The operator O4 then inserts energy, implying that the values thus obtained 
are unphysical. All these systematic errors can be studied in ChPT. Deviations due to the 
choice of unphysical masses and momenta already shows up at tree level [18], while all three 
systematic effects (including quenching) lead to one-loop contributions different from those 
calculated for physical masses and momenta in the unquenched theory [13]. In addition, 
at one loop one finds that there are power-like finite-volume corrections, which were also 
calculated in ref. [13]. 

The one-loop result for this unphysical matrix element in a finite volume (with peri- 
odic boundary conditions) and with mi = 1712 = my (hence = Mr, with M n the mass 
of the ud meson) for the partially quenched theory is, keeping only non-analytic corrections, 

{ n+A04\K + ) P h = Ml ff(l-N^ 2 lo g ^f (30) 
\ 1 41 lun P h V s V2/3 V (4vr/ P ) 2 6 A 2 K ' 



+ ^1 



14 



-31og-£ + F(M n L) 
Ap 



up to corrections vanishing faster than any power of L 1 , where the function F is given by 

^ , 17.827 12tt 2 

F (x) = + — • (31) 

X x° 

We note that the chiral logarithm due to the sea-quark loops does not diverge in the chiral 
limit my — > 0. Note also that the result Eq. (30) does not depend on the parameters itlq 
and a. For a discussion of contributions from 0(p 4 ) coefficients, see ref. [13]. 

The above result can also be derived from the "quark flow picture" [7,9], and the results 
of ref. [13]. First, take all quark masses equal, including the sea-quark mass. The quenched 
and unquenched results of ref. [13] then are the special cases obtained from Eq. (30) by 
setting N = and N = 3 respectively. The difference is due to the sea-quark loops which 
are present in the unquenched case, but not in the quenched case. Therefore, if we now 
have N instead of 3 sea quarks, one obtains the correct result by multiplying the difference 
between the unquenched and quenched results by N/3, yielding the first chiral logarithm in 
Eq. (30). But now we also identified which of the logarithms is due to sea-quark loops (the 
term proportional to N), and we conclude that we obtain the partially quenched result for 
mg 7^ my by replacing — > M^ s in the term linear in N. (Because of the structure of 
O4 and the 0(p 2 ) vertices from the effective lagrangian Eq. (3), there can be at most one 
sea-quark on the loop. Therefore, the result does not depend on M55 and the parameters 
mg and a, cf. Eq. (5).) Of course, one should keep in mind that the parameters f p and ctf 7 
depend on N. 

The power-like finite-volume corrections in Eq. (30) are independent of N. This follows 
from the fact that they originate from diagrams which do not contain sea-quark loops [13,21]. 

5. Numerical examples 



In this section, we will give some numerical examples of the differences between the 
one-loop estimates for and (7r + 7r°|04|i<' + ) calculated in the partially quenched theory 
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and in the "real world." For (7r+7r°|04|i ; r + ) we will also take the other systematic effects 
discussed in section 4 into account. We will take values for the lattice parameters typical of 
those used in recent numerical computations (for partially quenched results for Bg, see refs. 
[22,23]; for (7r + 7r°|04|.K' + ) we are not aware of any lattice data). 

The general strategy for our estimates is the same as in ref. [13]. We will set all 0(p 4 ) 
coefficients to zero. We choose the cutoffs A (for the full theory) and A p (for the partially 
quenched theory) to be 1 GeV or 770 MeV, independent of one another. The sensitivity 
under a change in A and A p is taken as an indication of the systematic error associated with 
our ignorance of the values of 0(p 4 ) coefficients. For the real- world values of f w , and 
m K we will use f n = 132 MeV, = 136 MeV and m K = 496 MeV. 

A. B K 

Bx in the full and partially quenched theories can be related by using Eq. (26) above 
and Eq. (36) of ref. [13]: 

R phys / f \ 4 

where 



P = ^, (33) 
1 + H 

with H = 0.724 (for A = 1 GeV) and H = 0.417 (for A = 770 MeV) is the numerical value 
of the relative one-loop correction for Bx in the real world, and H is the relative one-loop 
correction for the partially quenched theory in Eq. (26). To one-loop accuracy, M^ 2 can be 
replaced by M^. The factor P incorporates all one-loop corrections (P = 1 at tree level). 

Since the ratios f p / f and dnjc^i cannot be determined within ChPT, we will arbitrarily 
set them equal to one. This constitutes one of the major uncertainties of our method. 
Furthermore, we will take a = for simplicity. In Table 1, we list the numerical values of 
P for different combinations of cutoffs, choosing f p = f n , e = 1/2, Mgg = m^, N = 2, for 
various choices of the parameters y = {m^/?>)/M'g S and M^. 
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Ml 


y 


P (l) 
^(1) 


p (0.77) 
r (0.77) 


P (l) 
r (0.77) 


p (0.77) 
^(1) 


0.2 


0.5 
1 

oo 


1.02 
1.02 
1.01 


0.98 
0.98 
0.96 


1.19 
1.19 
1.17 


0.84 
0.84 
0.83 


0.4 


0.5 
1 

oo 


0.99 
0.98 
0.96 


1.13 
1.11 
1.05 


1.37 
1.35 
1.28 


0.81 
0.81 
0.79 


0.6 


0.5 
1 

oo 


1.11 
1.09 
1.03 


1.73 
1.64 
1.42 


2.10 
1.99 
1.73 


0.92 
0.90 
0.84 



Table 1. The factor P, Eq. (33), for different values of M^-, y and different combinations 
of A and A p . Other parameters are fixed at f p = f n , e = 1/2, M$g = mj( and N = 2. 
The superscript on P denotes A in GeV; the subscript on P denotes A p in GeV. is in 
GeV 2 . 

For M| = 0.6 GeV 2 , values for P as high as 2 are obtained, for Ap = 770 MeV. However, 
for this cutoff, the meson mass is probably too large for PQChPT to be reliable. Also, for 
= 0.6 GeV and y = 0.5, M^/m^, is of order one. We will therefore concentrate on 
the two lower masses in the following discussion. For these values, we see that P never 
differs from one by more than about 10%, if we choose A = A p , which is equivalent to the 
assumption that the bare parameters of the full and partially quenched theories are equal 
(this may not be unreasonable for N = 2). However, if we take the two cutoffs unequal, 
corrections can be as large as 20 — 30%. The results are fairly insensitive to changes in the 
value of y. 

B. K+ -> 7T + 7T 

The K + — > 7r + 7r^ matrix elements in the full and partially quenched theories can be 
related by using Eq. (30) above and Eqs. (43,87) of ref. [13] (replacing M\ x by M 2 and M\ s 
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by Mj> s ): 

(^V^+W = W J ^AOAKX nphv s ■ (34) 



with 

W 



(35) 



M2 



1 - log Tf + jfifr [-3 log ^ + F(M^)j 

where t/ = 0.0888 (for A = 1 GeV) and U = -0.0146 (for A = 770 MeV) is the numerical 
value of the relative one-loop correction in the real world. At tree level, W — 1. 

We list in Table 2a the numerical values of W for different combinations of the cutoffs 
A and A p , and for f p = f n , Mgg = m K , N = 2, several values of M% and a fixed volume 
L 3 such that M n L = 6 for M% = 0.2 GeV 2 . Here we may use M\ s = (M% + Af| s )/2. In 
Table 2b we list values of W for infinite volume, with all other parameters the same as in 
Table 2a. 

We see from Table 2 that one-loop corrections are always rather large, even for relatively 
small meson masses, and that sensitivity to the values of the cutoffs is significant. (For 
infinite volume, the corrections are not quite as big, cf. Table 2b.) This casts some doubt 
on the accuracy of one-loop ChPT in estimating the factor W, and one would expect that 
two-loop contributions are not small. As in the quenched case, the "correction factor" W is 
always substantially smaller than one. For a much more detailed discussion of uncertainties 
inherent to our estimates of such "correction factors," see ref. [13]. 

As can be seen from Eq. (30), the partially quenched result is closer to the unquenched 
case (for which N = 3, M\g = M\\) than to the quenched case. The large deviations from 
the tree-level value W = 1 are mostly due to the other systematic effects. For instance, in 
the unquenched theory, we find W = 0.56 (A = 1 GeV) and W = 0.57 (A = 770 MeV) for 



Ml = 0.2 GeV 2 . 
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0.2 0.59 0.60 0.66 0.53 

0.4 0.54 0.60 0.66 0.49 

0.6 0.55 0.66 0.73 0.49 



Table 2a. The factor W for different values of M% and different combinations of A and A p . 
Other parameters are fixed at: f p = f^, Mgg = m^, N = 2 and M n L = 6 for M% = 0.2. 
The superscript on W denotes A in GeV; the subscript on W denotes A p in GeV. M 2 is in 
GeV 2 . 



M 2 


< 


^(0-77) 
vv (0.77) 


W (0.77) 




0.2 


0.68 


0.71 


0.78 


0.62 


0.4 


0.65 


0.75 


0.83 


0.59 


0.6 


0.68 


0.90 


0.99 


0.62 



Table 2b. The factor W for different values of Ml and different combinations of A and A p . 
Other parameters are fixed at: f p — f n , Mgg = m^, N = 2 and M n L = oo. The superscript 
on W denotes A in GeV; the subscript on W denotes A p in GeV. M 2 is in GeV 2 . 

6. Conclusion 

In this paper, we applied PQChPT to the calculation of the quark-mass dependence 
of various quantities to one loop, restricting ourselves to degenerate sea-quark masses. For 
Goldstone boson masses and decay constants, we extended results obtained earlier in ref. 
[4]: we investigated the sensitivity of the chiral logarithms to the singlet part of the rj' mass, 
niQ. Since the "remnants" of the rj' are an essential part of the fully quenched theory, this 
is a natural question to ask in the partially quenched theory. We found that the coefficients 
of the chiral logarithms for the Goldstone meson masses are sensitive to the value of mg for 
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typical values of the sea-quark mass. The decay constants depend on mg to a lesser extent, 
which is related to the fact that in the limit of degenerate quark masses, the axial current 
does not couple to the rj . 

We also calculated one-loop contributions to the chiral condensate, Br and the K + — > 
7r + 7r decay amplitude, extending earlier work in QChPT [9,16,13] to the partially quenched 
case. For K + — > 7r + 7r° we only considered the case of degenerate valence-quark masses. 

We considered some numerical examples of the comparison between and the K + — > 
7T + 7r° decay amplitude in the real world and in partially quenched QCD with values of the 
parameters typical of current lattice computations. For B^, we found that, for small enough 
meson masses, the "correction factor" may be very close to one, but with an uncertainty 
which could be as large as 20 — 30%. For K + — > 7T+7T , we also took into account that 
typical lattice computations are done at unphysical (degenerate) valence-quark masses and 
external momenta, and we included the leading finite-volume corrections. We found that 
the "correction factor" is always much smaller than one (of order one-half, with large uncer- 
tainties). This is mostly due to the unphysical choice of masses and momenta, and not to 
partial quenching. 
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